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Preface 


This book introduces nonlinear control systems for control engi- 
neering and science to graduate, undergraduate students and re- 
searchers; it targets control engineering students who do not like 
to do not have time to derive and prove mathematical results for 
nonlinear control systems. It can be serve as a text book for nonlin- 
ear control systems, especially for feedback linearization techniques 
which is a common approach in controlling nonlinear systems. 

The development of computer software for nonlinear control 
systems has provided many benefits for teaching, research, and the 
development of control systems design. MATLAB® is considered 
the dominant software platforms for linear and nonlinear control 
systems analysis. This book contains a MATLAB-based program 
that helps teachers, students, and researchers for finding feedback 
linearization controller for a specific nonlinear class of systems. 

Chapter 1 is an introduction to theoretical aspects of nonlinear 
feedback linearization techniques. We use simple and complex ex- 
amples to better illustrate the method on how to find the feedback 
linearization controller for single and multi input output systems. 
In section 1.9, a MATLAB-based program is developed and used 
to find feedback linearization control using symbolic MATLAB li- 
braries for a special class of nonlinear systems. 

Chapter 2, the concept of structure variable control will be in- 
troduced in intuitive way, with illustrative examples that makes the 
reader familiar with the theory of sliding mode control in an easy 
and simple way. After having defined the control law, the param- 
eters of the controller should be identified, hence a simple method 
is presented to find those parameters. In section 2.4 a MATLAB- 
based function is developed to find the sliding mode surfaces and 
controllers using symbolic MATLAB library for SISO and MIMO 
nonlinear dynamical systems that could be written in a special 
form. 


vii 
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This textbook provides an easy way to learn nonlinear control 
systems such as feedback linearization technique and sliding mode 
control (structrue variable control) which are the most used tech- 
niques in nonlinear control dynamical systems; therefore, teachers, 
students and researchers are all in need of learning how to handle 
such techniques since they are too difficult for them to handle such 
nonlinear controllers, especially for more complicated systems such 
as induction motor, satellite, and vehicles dynamical models. Thus, 
this document is an excellent resource for learning the principle of 
feedback linearization and sliding mode techniques in an easy and 
simple way. This book: 


(a) provides a briefs description of the feedback linearization and 
sliding mode control strategies 


(b) includes a simple method on how to determine the right and 
appropriate controller (P-PI-PID) for feedback linearization 
control strategy 


(c) provides a MATLAB symbolic based programs that can solve 
the Lie derivative of any nonlinear dynamical system that 
could be written in a special form 


(d) provides a symbolic MATLAB-based function for finding the 
feedback linearization, and sliding mode controllers are devel- 
oped and tested using several examples 


(e) introduces a simple method for finding the approximate sliding 
mode controller parameters 


The relative degree of the system which can be found using the 
provided MATLAB code can be used to determine the degree of 
the Sliding Mode Surface, which is a very important scalar for 
those who are working with Variable Structure Control. 

There are many examples used in the last chapter with a mod- 
ified MATLAB program such as: 


(a) A nonlinear dynamical model of a pendulum system 
(b) A nonlinear dynamical model of Ven der pol system 
(c) A nonlinear dynamical mode of an aircraft 

) 


(d) A nonlinear dynamical model of an induction motor 


Preface Mix 


(e) A nonlinear dynamical model of a Permanent magnet syn- 
chronous motor. 


Where the program used to construct the nonlinear controller 
uses symbolic computations such that the user should provide the 
program with the necessary functions f(x), g(#) and h(a) using 
the symbolic library. 

Note: Make sure the function syms exists on your MATLAB 
by typing help syms on MATLAB. 

MATLAB® is a trademark of The MathWorks, Inc. and is used 
with permission. The MathWorks does not warrant the accuracy 
of the text or exercises in this book. This book’s use or discussion 
of MATLAB® software or related products does not constitute 
endorsement or sponsorship by The MathWorks of a particular 
pedagogical approach or particular use of the MATLAB® software. 
MATLAB® is a registered trademark of The MathWorks, Inc. For 
product information, please contact: 

The MathWorks, Inc. 

3 Apple Hill Drive 

Natick, MA, 01760-2098 USA 
Tel: 508-647-7000 

Fax: 508-647-7001 

E-mail: info@mathworks.com 
Web: www.mathworks.com 


Taylor & Francis 
Taylor & Francis Group 


http://taylorandfrancis.com 


CHAPTER 1 


Feedback 
Linearization Control 


FEEDBACK LINEARIZATION is a common approach used in 

controlling nonlinear systems. The approach involves coming 
up with a transformation to the nonlinear system into equal linear 
system that could be controlled easily using a new input control. 
Feedback linearization could be applied to nonlinear systems of the 
form: 


& =f(x) + g(a)u (1.1) 
y =h(2) 


Where x € #” is the state vector, x € R? is the vector of inputs, 
and y € &™” is the vector of outputs. The goal is to develop a 
control input: 


u = D(x)~" [—A(z) + 9] (1.2) 


That renders a linear input-output map between the new input v 
and the output. An outer-loop control strategy for the resulting 
linear control system can then be applied. 


1.1. FEEDBACK LINEARIZATION OF SISO SYSTEMS 


Let’s consider the following SISO nonlinear dynamical system de- 
scribed by: 


&= f(x) + gu 
y = h(z) (1.3) 
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Where x = [21, £2, ...., Un] is the state vecotr; u input vector to the 
system; f(a) and g(x) are infinite differentiable vectors; hence the 
input-output feedback control law u is defined for SISO nonlinear 
systems in (1.6) by the following relation: 


1 


u= Hina +v) (1.4) 


Where r is the relative degree of the system; L,L"~' f(x) and 
L'-h(x) are the lie derivative of h(x) along the vector fields g(x) 
and f(a), respectively. The notion of the lie derivatives and how 
its calculated will be given in the subsequent sections. Thus, the 
overall system control law u and nonlinear system (1.6) combined 
in series yields to a linear system with order r (relative degree) 
from the new input v to output y 


y’ =v (1.5) 
Where v is the new input vector to the system; that will be used 


as new controlled input. 


RELATIVE DEGREE OF SISO SYSTEMS 


1. r=n the system admits an exact state linearization 
2. r <n the system admits partial feedback linearization 


3. r > n the system does not admit an input output feedback 
linearization 


1.2 FEEDBACK LINEARIZATION OF MIMO SYSTEMS 


Feedback linearization control for MIMO system is applied to dy- 
namical systems that have the following form: 


&=f(x)+ ye Gi(X) Us 


y1 =hy (2x) 
y2 =ho(x) (1.6) 


Feedback Linearization of MIMO Systems @ 3 


Equation 1.6 is called a square system, where the number of inputs 
is equal to the number of outputs. 

The control law for MIMO systems can be obtained using lie 
derivative. 


yi 
¥2 | = A(x) + D(x)u (1.8) 
vr 
Where hy(2) 
Lh; (a 
pred 
A(x) = oe (1.9) 
Lie” hm (2) 
Dg le). Lg ihe) 
D(a) = (1.10) 


Lig tle) Lgl ea) 


Where D(z) is called the decoupling matrix of the system, and it 
is not always nonsingular, and therefore the linearization control 
is then achieved if and only if the decoupling matrix is invertible 


L''hj(x) The lie derivative of h(x) along the vector field f(z) 
Lom Le hace) The lie derivative of h(a) along the vector field 
g(x) 


Therefore, the linearization control law is then written in com- 
pact form as: 
u = D(x)~" [—A(z) + 9] (1.11) 


Where v is the new input control vector which has the same di- 
mension as u. 


RELATIVE DEGREE OF MIMO SYSTEMS 


lo r=rjt+ret....1m =n then the system admits an exact 


2.r=7r+rgt..%m < n then the system admits partial 
feedback linearization 


3. r > n the system does not admit an input output feedback 
linearization for both types of systems 
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1.3 RELATIVE DEGREE 


The relative degree of a dynamic systems is obtained by differen- 
tiating the output y until the input u appears 


Theorem 1 A nonlinear system is said to have a relative degree 
r at x off: 


LL *h(a) #0 for r#An and r<n 


Where Lh is the lie derivative of h along the vector field f(x) 


Where Lrh is the derivative of h along the vector field g(x) 
i=1 


gh(a) = SO ge) 


n 


1.4 MATLAB® PROGRAM DESCRIPTION 


The technique of an input-output feedback linearization is used for 
many control applications, which transforms a complicated and 
coupled nonlinear system into an equivalent linear system that 
could be controlled using linear control system tools, such as pole 
placement technique. Thus, to find such control law a MATLAB 
script is developed to facilitate the tedious computations and avoid- 
ing errors due to a complexity of the system. 


1.5 _PROBLEM FORMULATION 


Let’s consider the following nonlinear mathematical system repre- 
sentation: 


= fle) t+ Or? gi(z)ui 
ae 2 


Where x € &” is the state vector; x € R” is the input vector of 
the system; h;(x)a € ¥ Is an analytical function of x; f(x) and 
g(x) are infinitely differentiable vector fields. 
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1.5.1. Programme flow chart 


The flow chart of Figure 1 gives a detailed explanation of how the 
program works; therefore, the following steps are used to obtain 
the nonlinear control law using the developed program: 


(a) Introduce the nonlinear functions f(x),g(x) and h(a) as a sym- 
bolic functions to the program 


(b) Find the lie derivative along the vector fields f(a) and g(x) 


(c) If the lie derivative along the vector fields g(x) = 0 and the 
number of differentiation is less than n (the system state num- 
ber), then set f(a) = Lrh(x) and repeat step 2 


(d) If the lie derivative along the vector fields Lgh(x) # 0 and the 
number of differentiation is or equal to n (the system state 
number), then prints the lie derivative along the vector fields 
f(x) and g(x) as well as the differentiation number r, which is 
the relative degree, which could be used to identify the con- 
troller type 


— 
fo) 
——- 


If the lie derivative along the vector fields L,h(x) = 0 and the 
number of differentiation is equal or greater than n (the system 
state number), then the system does not admit an input output 
feedback linearization 


(Start ») 


{| 1.2). he) || 
<“alculate the Lie ~ 
/ derivatives Lhf.c / 


Yes - 
Lhe—0 — fix<)—Lhbhft Z 


No + 
. . 
) Output LbfL bs| ; 
J 


4 Stop 


Figure 1.1 The flow chart of the proposed algorithm 
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1.6 


MATLAB CODE OF LIE DERIVATIVE 


The aim of feedback linearization is to produce a transformation 
system whose states are the output y and its derivatives. Therefore, 


the 


lie derivative is then used, and since it’s an important tool for 


analysis and synthesis, a MATLAB program has been developed 


for 


any class of nonlinear system that could be written in the form 


of an equation (1.12). The two programs work as follows: 


I 


II 


— 


oe 


ae 


iG 
h 
x 
a 


oP ol? 


oe 


The LieDerivative(h, x) program finds the partial derivative 
of h(a) along x; as described by the following equation: f(a) = 
On a) jy. thus, the user should provide the program with the 
vector field 


The solvelieder(Lh, fx, g) gives the lied derivatives of the vec- 
tor fields f(z) and h(x) along the vector fields f(x) and g(x) 
respectively. Where the user has to provide the program with 
the function f(a) and g(a) as well as the output derivative 
from the function LieDerivative(h, x) 


ction df=LieDerivative (h, x) 


The LieDerivaive MATLAB function is used 


o find the jacobian vector of a given output 
(32) : Is the output function 

: The state vector 
£ : The jacobian of h along x 


if nargin<2 & nargin== 


end 
df= 


error('not enough input argument"); 


[li 


n=length (x) ; 


for 


end 
af; 
end 


fun 


AP JP oP oP ol? 


ae 


ii=lin 

XX=x (ii); 
dff=diff (H, xx); 
df=[df,dff]; 


ction [lhf lhg]=solvelieder (Lh, £x,g) 

The solvelieder MATLAB function is used to find 
the lie derivatives of the functions f(x) and g(x) 
along the vector field h(x) 

Lh : The jacobian vector of h along x 

tx : The function f(x) 

g : The input function g(x) 
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LHg=[1]; 

Lhf=Lh+fx; 

[n,b]=size(G); 

for ii=l:b 
Lgh=Lh+*G(:,ii); 
LHg=[LHg, Lgh]; 

end 

lhg=LHg; 


1.6.1 Example for SISO systems 
Find the lie derivative of the following system with h(x) = x1 
{ pa Nie (1.13) 


tg = —sin(a1)+u 


Hand Calculation The lie derivative of any nonlinear dynamical 
systems that has the form of Equation 1.1 is obtained by differen- 
tiating the output of the system; hence, the first lie derivative of 
Equation 1.13 is given by: 
. Oh(a) 
oe 
=[1 0] (F(@)+9u) (1.14) 


Therefore, from Equation 1.14 the first lie derivative is given by: 


p40] ( eee ) (isl) 
y =Lyh(z) = x2 (1.16) 


£ 


MATLAB Computation The first lie derivative of equation (1.13) 
is then obtained using MATLAB program described on the previ- 
ous section as follows: 


syms xl x2 u % symbolic presentation 
fx=[x2;-—sin(x1)]; 

g=[0;1]; 

x=[xl x2];H=x1; 
LH=Lie_Derivative (H, x) 

[lhf lhg]=solve_lie_der (LH, fx,g) 


Program output 
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LH = [ 1, 0] 
lhf = x2 
lhg =0 


The first lie derivative of the dynamical system in Equation 1.13 
is displayed using MATLAB as Lrh(x) = Ihf and Lgh = lhg. 


1.7 CONTROLLER TYPES 


After having defined the matrix A(a#) and D(x) as well as the rel- 
ative degree of the system, which has a significant importance of 
finding the appropriate type of controller, and since that the trans- 
fer function of the linearized system is given in terms of the relative 
degree are as: 


y il 


U s” 


(1.17) 


Hence, the controller type will be defined in an appropriate manner 
using values of the relative degree are as: 


e If r=1 then the appropriate controller is a Proportional P 


e Ifr=2 then the appropriate controller is Proportional Deriva- 
tive PD 


e If r=3 or greater than that, the appropriate controller is 
Proportional-Integra-Derivative PID 


1.8 FEEDBACK LINEARIZATION CONTROLLER 
EXAMPLES 


An illustrative examples will be used to show users how to use 
MATLAB code for finding the nonlinear control law. Therefore, 
for better readability of the program, nonlinear dynamical systems 
are used. 


1.8.1 Example for SISO systems 


The nonlinear dynamical model of a pendulum is written in com- 
pact form as 


é= f(x) +gu= ( ies) ) + | : J (1.18) 


Find the control law u that linearized model with h(x) = x2? 
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1.8.1.1. Solution Using Hand 


Following the steps mentioned in the flow chart of Figure 1, then 
the control law will be obtained by differentiating the output h till 
u appears as: 


_ h(x), _ Ah(x) 


u(z) = =  (F (2) + ou) (1.19) 

Where 
Lth(z) = OME) F() =(0 1] | enced =~—sin(x,) (1.20) 
DA(a) = On) a(x) =1(0 1) | : | ay (1.21) 


Where the superscript on Equation 1.17 denotes the relative degree 
of the system which is 1, and Lh, Lgh are the lie derivatives of 
the vector field h along the vector fields f(x) and g(x), respectively. 
Furthermore, the control law u is obtained using Equation 1.9 as: 


u = sin(a1) +0 (1.22) 


1.8.1.2 Solution Using MATLAB 


Equation 1.22 is then obtained using MATLAB script presented at 
the end of this document by: 


The Nonlinear systems should be written in the following 

form—Feedback Linearization Controller for a class of 
Nonlinear systems—State space equations 
x=f(x)+g(x)u 


Parmeters 
parameters = 


Empty string: 1—by—0 


Number of states:=2 

Number of inputs:=1 

The vector f(x) :=[x2;—sin(x1) ] 
The vector g(x) :=[0;1] 

The output variables:=x2 

The relative degree of hl 
equal:=1 
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The feedbacklinearization controller Uc: 
Uc = 
ul + sin(x1) 


By comparing the control law obtained in Equation 1.22 using 
hand calculation with that obtained using MATLAB script code 
shows the effectiveness of the program for single input-single out- 
put system. 


1.8.2 Example for MIMO systems 


In order to show the effectiveness of the program for a system that 
has a multi-input-multi-output, the following nonlinear dynamcial 
system will be used: 


Vy + L122 4 1 0 
—sin(x1) 0 1 


“&= f(x)+gu= | 
ae hy 
The output vector of the system is given by: h(a) = = 
XY 
X2 
1.8.2.1 Soluting Using Hand 


Following the steps mentioned in the flow chart of Figure 1; then 
the control law will be obtained by differentiating the output h till 
u appears as: 


_ h(x), _ Ah(z) 


ia) = Os = BO He) 49x) (1.24) 

Where 
h(a) = On) F(a) ie) | nies = 2, +a (1.25) 
3ho(2) = ae) fla) =[0 1] | Sua = —sin(a2) (1.26) 
L1hy (x) = OM) (0) ={1 0] | ; : es (1.27) 
ryta(a)= "Fra =o a] p]=2 38) 
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Where the superscript on Equation 1.10 denotes the relative degree 
of the system which is 1, and Lyh, Lgh are the lie derivatives of 
the vector field h along the vector fields f(x)and g(x), respectively. 
Furthermore, the control law u is obtained using Equation 1.11 as: 


Pe ]=[5 3][egte] om 


1.8.2.2 Solution Using MATLAB 


Equation 1.29 is then obtained using the MATLAB script pre- 
sented at the end of that document as shown: 


The Nonlinear systems should be written in the following 

form—Feedback Linearization Controller for a class of 
Nonlinear systems—State space equations 
x=f(x)t+g(x)u 


Parmeters 
parameters = 
Empty string: 1—by—0 


Number of states:=2 

Number of inputs:=2 

The vector f(x) :=[x1+x1*x2;—sin (x1) ] 

The vector g(x):=[1 0;0 1]; 

The output variables:=[x1; x2] 

The feedback linearization controller Uc: 
inv (Lhg) + (-Lhf+u) 

The relative degree of hl 


equal:=1 

The relative degree of h2 
equal:=1 

Lhf = 


xl + xl*x2 


—sin (x1) 
Lhg = 
[ 1, 0] 
[ 0, 1] 


By comparing the control law obtained in Equation 1.29 using 
hand calculation with that obtained using MATLAB script code 
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shows the effectiveness of the program for multi-input-multi-output 
systems. 


1.9 MATLAB BASED FUNCTION: FEEDBACK 
LINEARIZATION 


EEDBACK LINEARIZATION problems are then solved using 

MATLAB-based function; and for better readability of the 
programmes illustrated on the previous chapters, we will write 
them on another manner that simplify the use for users; hence, 
the program that should be used by the user will be written as: 


1.10 MATLAB TEST FUNCTION 


clear all;clce 
disp (' ')i 
disp('The Nonlinear systems should "he 
disp('be written in the following form Pe 
disp('State space equations x=f (x) +g(x)u'); 
disp (' ')i 
% The your system contains 
% Input the extra parameters 
par=input ("Parameters ','s'); 
eval(sprintf('syms %s',par)); 
parameters=sprintf('%s',par) 
%% Declare how many states and inputs 
n=input('Number of states:="'); 
nin=input ('Number of inputs:='); 
x=sym(zeros(1,n)); 
u=sym(zeros(1,nin)); 
for j=lin 
eval(sprintf('syms x%d',4)) 
x(:,j)=sprintf('x%d',j); 
end 
for k=1l:nin 
eval (sprintf('syms u%d',k)); 
u(:,k)=sprintf ('u%sd',k); 


end 

% Enter the functions from the keyboard 
f=input ('The vector f(x):=','s"'); 
g=input ('The vector g(x):=','s'); 
He=input ('The output variables:=','s'); 


Represent all the functions 

S£(x), g(x) and h(x) on a symbolic format 
fx=sym(f); 

g=sym(g) ; 

Hc=sym (Hc) ; % 
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% Use the inoutfeedbacklinearization.m 
% programm to generate the desired functions 
[Lhf Lhg]=inoutfeedbacklinearization(fx,g,Hc,x) 


“newlabell.101537 
1.10.1 Feedback linearization MATLAB function 


Thus the inoutfeedbacklinearization.m function is then written 
using a nested MATLAB functions instead of using the programs 
separately as shown on the previous chapters; therefore the whole 
code is given by: 


function [Lhf Lhg]=inoutfeedbacklinearization(fx,g,h,x) 

The function inoutfeedbacklinearization is used to find the 
the feedbacklinearization control law for SISO and MIMO 
nonlinear systems using symbolic MATLAB library; 

The user should provide the program with the following 


% inputs 

ei : The system function f(x) 

Sg : The system output function g(x) 

%& h(x) The vector of outputs h(x)=[x1l;x2, He-4 ull 

ane 4 : The state vector x=[x1,x2, ,xn] 

% After having provided the program the necessary input 

% functions 

% The program will output the following variables 

% Lhf : The lie derivative of h(x) along the function f (x) 
% Lhg : The lie derivative of h(x) along the function g(x) 
% which is called the decoupling matrix 

% u : The vector of inputs u=[ul;u2;...;un] 

% The control law will be given by the following formula 

% u= inv (Lhg) «(—Lhgtv) 


if nargin <4 
error('Not enough input argument'); 
end 
k=1; 
Lhg=[];Lhf=[]; 
nb=length (h) ; 
while k<length(h)+1 
hi=h (k) ; 
for i=l:nbt1l 
% this Lie derivative function 
df=Lie_Derivative(hl,x); 
% solve for the g 
[lhf lhg]=solve_lie_der(df,fx,g); 
[n b]=size(lhg); 
for ii=l:n 
d=any (lhg(ii,:)~=0); 
end 
if d==1; 
disp(['The relative degree of h',num2str(k)]), 
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disp(['equal:=',num2str(i)]); 
break; 
else 
hl=lhf; 
end 


if i==nb+1 ¢& d== 

disp(['The system dose not admit NFL']); 
return; 

end 
end 
Lhg=[Lhg; lhg]; 

Lhf=[Lhf;lhf]; 

k=k+1; 

end 

function df=LieDerivative (H, x) 

The LieDerivative MATLAB function is used 

to find the jacobian vector of a given output 
H (x) : Is the output function 

The state vector 

The jacobian of h along x 


* 


AP PoP ol? lO 


©: 
Fh 


if nargin<2 & nargin== 
error('not enough input argument'); 
end 
df=[]; 
n=length (x) ; 
for ii=l:n 
XX=xX (ii); 
dff=diff (H, xx); 
df=[df,dff]; 
end 
df; 
end 
function [lhf lhg]=solvelieder (LH, fx,G) 
The solvelieder MATLAB function is used to find 
the lie derivatives of the functions f(x) and g(x) 
along the vector field h(x) 


AP © oP oP a OP 


LH : The jacobian vector of h along x 
fx : The function f(x) 
g : The input function g(x) 

LHg=[]; 

1lhf=LH+fx; 

[n,b]=size(G); 

for ii=l:b 


Lgoh=LH*G(:,ii); 
LHg=[LHg, Lgh]; 
end 
lhg=LHg; 
end 
end 
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In order to use the two programmes, the user should follow the 
steps: 

e Save the function inout feedbacklinearization.m given on 


the code above as shown on Figure 1.2 


DSe|4eOC|SB-|Medhe|B- HH BOM DB) stace| Bae 


GSi1 function [Lhf Lhg] =inoutfeedbacklinearization(£x,g,h,x) 


46 % the vector field func(x) 

47 &H 3 the output vector 
46 tx : the state varibales 
a9 


if nargin<2 & nargin==0 
error('not enough input argument"); 

end 

af=[]- 

n=length (x) ; 

for ii=1:n 


w 
o 
| UO TRF ae tea Ta SR eet Ci Sic aR | 


xxX=x (ii)? 
57 dff=ditt(H,xx); 
58 af=[df,dff];: 
59 end 
60 af: 
61 end 
62 function [1hf 1hg]=solvelieder (LH, fx,G) 
63 % This equation solves the Lie derivitavies that is mulitiplied by 
64 % The output vectors Lhg=Lh*g*u where 
65 & LH : The Lie derivitavies of h(x) along the vector field fx 
66 % fx : The vector field f(x) that describe the system 
67 &G : The vector filed of the input 
68 — LHo=[]- 
69 - Lnf=LH*fx; 
70 - [n,b]=size(G); 
7 = for ii=1:b 
72- Lgh=LH*G(:,ii); 


Figure 1.2 inoutfeedbacklinearization MATLAB Code 


e The call program where the user has to provide the program 
with necessary input functions is given as 


clear all;clce 

disp (' ')i 

disp('The Nonlinear systems should be written in the 
following form'); 


disp ('—Feedback Linearization Controller for a class of 
Nonlinear systems—'); 

disp(' State space equations x=f (x)+g(x)u 

")i 

disp (' ')i 


20 


%% Declare how many states and inputs 
% The your system contains 

% Inout the extra parameters 
par=input ('Parameters ','s'); 
eval(sprintf('syms %s',par)); 
parameters=sprintf('%s',par) 

n=input ('Number of states:='); 
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nin=input ('Number of inputs:='); 
x=sym(zeros(1,n)); 
u=sym(zeros(1,nin)); 
for j=lin 
eval(sprintf('syms x%d',4)) 
x(:,j)=sprintf('x%d',j4); 
end 
for k=1l:nin 
eval (sprintf('syms u%d',k)); 
u(:,k)=sprintf ('u%d',k); 


end 

6% Enter the functions from the keyboard 
f=input ('The vector f(x):=','s'); 
g=input ('The vector g(x):='); 

He=input ('The output variables:=','s'); 


%% Represent all the functions f(x), g(x) and h(x) ona 

6% symbolic format 

fx=sym(f); 

sg=sym (gq); 

Hc=sym (Hc) ; % 

%% Use the inoutfeedbacklinearization.m programm to generate 
%% the desired functions 

disp(['The feedbacklinearization controller Uc:']); 
disp(['inv(Lhg) *(—Lhf+u) ']) 

[Lhf Lhg]=in_out_feedback_linearization (fx,g,Hc,x) 


1.11 ILLUSTRATIVE EXAMPLES 


1.11.1 Aircraft altitude dynamics 


Let us determine the associated internal dynamics of the aircraft 
by defining the state x as [a, wh, h], thus, the equations of motion 
can be written as: 


Ly =22 

to = — 40g — 47, 4+ 3E 

t3 =@4 (1.30) 
t4 =62, —E (1.31) 


Where E is the input vector, a = 2, and a = £2 are the angle and 
angular speed of the aircraft, respectively; h is the altitude and h is 
the altitude speed see Figure 1.3. By applying the MATLAB-based 
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Figure 1.3 Dynamic charactercistics of an aircraft 


program the feedback linearization controller of the aircraft is given 
below: 


The Nonlinear systems should be written in the following 

form—Feedback Linearization Controller for a class of 
Nonlinear systems—State space equations 
x=f(x)+g(x)u 


parameters = 
Empty string: 1—by—0 


Number of states:=4 

Number of inputs:=1 

The vector f(x) :=[x2;—4*x2—4«x1; x4; 6*x1] 
he vector g(x) :=[0;3;0;—-1]; 

he output variables:=x3 

he feedbacklinearization controller Uc: 
inv (Lhg) * (-Lhf+u) 

The relative degree of hl 

equal:=2 


Lhf = 


From the output of the program shown above and Equation 1.9 
the new control vector v is given by: 


E=62;-—v (1.32) 
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1.11.2 Asynchronous motor speed control 


The dynamical model of an induction machine is given by: 


— ye. + #3 + pKar5x4 


—Y@2 47 704 — pK25%x3 
f(z) = Atty = ts — Pets (1.33) 
rx = Ta + pr5x3 


Por, (wate — tae) — Fears — Fh 


1 T 
0 0 0 
oL 
g= 7 1 (1.34) 
0. 2. U4 
: = M? 
Where the parameters o,7,7; are defined as) of = 1 —- Tr 


2 rf . . 
Y= tee feat T, = ie Where o is the scattering coefficient, 


T,. is the time constant of the rotor dynamics, J;, is the rotor in- 
ertia, fm is the mechanical viscous damping, p is the number of 
pole pairs, Tq is the unknown external load torque which will be 
estimated as well. 

The state variables isa, isq, Grd, Usd, Usq, are the stator currents, ro- 
tor flux linkages, stator terminal voltages respectively. and for the 
parameters L,, L,, M, R, are the rotor inductance, stator induc- 
tance, mutual inductance, stator resistance, and rotor resistance, 
respectively 


U= [Use Use i [ese 158 Gra orp Q] (1.35) 
The output vector is then given by: 


h(a) = | ae | (1.36) 


Using the program, the output will be: 


The Nonlinear systems should be written in the following 

form—Feedback Linearization Controller for a class of 
Nonlinear systems—State space equations 
x=f(x)+g(x)u 


Parameters Tl gamma K Tr p fm Jm Lr M Ls sigma 
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parameters = 
Tl gamma K Tr p fm Jm Lr M Ls sigma 


Number of states:=5 
Number of inputs:=2 
The vector f(x) :=[—gamma*x1+ (K/Tr) *x3+p*K*x5*x4; 
—gamma*x2+ (K/Tr) *x4—p*K*x5*x3; 
(M/Tr) «x1—(1/Tr) *x3—p*x5*x4; (M/Tr) «x2 
—(1/Tr) *«x4+p*«x5«x3; 
px (M/ (Jm«Lr) ) * (K3*x2—x4*x1) 
—(fm/Jm) *x5—T1/Jm] 
The vector g(x):=[1/(sigmaxLs) 0;0 1/(sigmaxLs);0 0;0 0;0 0] 
The output variables:=[x3°2+x4°2;x5] 
The feedbacklinearization controller Uc: 
inv (Lhg) * (-Lhf+u) 
The relative degree of hl 


equal:=2 

The relative degree of h2 

equal:=2 

Lhf = 

((4*x3)/Tr — (2*M*x1)/Tr) *(x3/Tr + p*x4*x5 — (M*x1)/Tr) 
— ((4*x4)/Tr — (2*M*x2)/Tr)* (p*x3*x5 — x4/Tr + (M*x2)/Tr) 
+ (2*M*x3* ((K*x3)/Tr — xl*gamma + Kxp*x4«x5))/Tr 

— (2*M«x4« (x2*gamma — (K*x4)/Tr + K*xpx*x3*x5))/Tr ; 


(fmx« (T1/dm + (fm*x5)/JIm + (M*px (xl*x4 — x2*x3))/(Im*Lr)))/dm 
— (M*p*x3* (x2*gamma — (K*x4)/Tr + K*px*x3*x5)) / (JIm*Lr) 
— (Mxp*xx4« ((K*x3)/Tr — xlxgamma + K*xpx*x4«x5)) / (Jm*Lr) 
— (Mep*xxlx* (p*x3%*x5 — x4/Tr + (M*x2)/Tr))/ (Jm*Lr) 
(Mxp*x2* (x3/Tr + pxx4*x5 — (M*x1)/Tr))/(Jm«Lr) 


Lhg = 


[ (2*M*x3) / (Ls*Tr*sigma), (2*M*x4) / (Ls*Tr*xsigma) ] 
[ —(M*p*x4) /(dm*xLr«Ls*sigma), (M*xp*x3) /(Jm*Lr*Ls*sigma) ] 


Taylor & Francis 
Taylor & Francis Group 


http://taylorandfrancis.com 


CHAPTER O 


Sliding Mode Control 


N practical control problems there will be a mismatch between 

the real plant and its developed mathematical model used for 
control purposes. These mismatches arise from uncertainties due 
to changes in the plant parameters or due to external disturbances. 

The performance of any nonlinear dynamical systems that have 
such uncertainties becomes a challenging task for control engineer- 
ing; a robust control method is developed to solve this problem 
which so-called Sliding mode control or Structure variable 
control [8] is used to construct the sliding mode surface; and, 
therefore the sliding mode controller for a class of nonlinear dy- 
namical systems that can be written as: 


2.1 SLIDING MODE CONTROL THEORY 


Let’s consider the following nonlinear dynamical system:. 


& = f(z,t)+9(a,t)u (2.1) 
y =h(a,t) (2.2) 


Where y and u denote the output and input variable, x is the 
state vector. The sliding mode control makes the output y tracks 
a desired input reference y,, so that the system motion is kept on 
the manifold $, which is defined by: 


S =o(a,t) =0 (2.3) 
S represents the sliding mode surface; and, therefore, the state 


variables has to be brought to the siding surface in order to obtain 
the control law u; which is given by: 
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Ueq : is the equivalent control vector, which is obtained from the 
derivative of Equation 2.3 


the correction factor) vector which 


S'), Where k, is a controlled gain, 
= 


Un : is the discontinues control ( 
is given by: un = —kpsgn( 
sgn is the function sign(S 


The sliding mode surface S form takes the following form: 
S=o=(—+k)" Te (2.4) 


Where r is the degree relative of the system, e is the errors between 
the output vector and the desired reference input; k is a positive pa- 
rameter which can be chosen arbitrarily or using a simple method 
that could lead to the appropriate choice. 


2.1.1. SISO Sliding Mode Control 


The sliding mode control law for SISO systems can be obtained 
using the equivalent and switching control law as follows: 

The equivalent control law: is given by setting the derivative 
of the sliding mode surface to zero as follows: 


S$=0 (2.5) 


And then solve for the control input law u which will be the equiv- 
alent control 
The switching Control is given by following equation: 


Un = S = —kpsgn(S) (2.6) 


Therefore, the overall control law u = Ueg + Un is then given by 
setting: 

dd wa 

TF + bye) = -kpsgn(S) (2.7) 


And solve for the control law u, the idea is better illustrated using 
examples 
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2.1.2 MIMO Sliding Mode Control 
Consider the MIMO nonlinear system given by: 


b= f(a) +0 gui 
w=1 


yi =hi(x) 
Yy2 =h2(z) (2.8) 
The nonlinear system (2.8) has a vector relative degree 14, ..... Tn 


and has a vector output h;. We assign a sliding mode surface for 
each output; a sliding mode control for MIMO systems are obtained 
by following the steps outlined for SISO systems. 

The MIMO sliding mode control law is obtained by setting 
the derivatives of the sliding mode surfaces equal to the switching 
surfaces as shown: 


$1 = — Kp,sgn(S1) 
Sy = — Ky, sgn(S2) 

(2.10) 
Sn=- Kp, 89n(Sn) 


substituting the sliding mode surface given in (2.4) onto (2.10) 
gives: 


“(e + ky)"~*e1) = — Kp, sgn(51) 
“(4 + kz)" *e2) = — Kp, sgn(S2) 

(2.11) 
“(e sitet) cise sets) 


Where ky, ko,...kn and ky, kp2,...kpn are a controlled parameters, 
€1, €2,---Cn are the output errors of the system. The control law is 
then obtained by solving Equation 2.11 for the vector control u;. 
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2.2 SLIDING MODE CONTROL EXAMPLES 
2.2.1 Vander pol system 


Consider the Van der Pol dynamical equations [4]: 


Ly =% 
bo = —a, +e(1—2?)ro +u (2.12) 
Yr21 


Find the sliding mode controller for the Van der Pol system. 


2.2.1.1. Solution Using Hand 
The sliding mode surface is then defined by: 


+k)"‘e (2:13) 


Where the error e = (#1, — 21) is written in terms of the input 
reference x, and the output state x1; k is a controlled parameter; 
the relative degree r = 2 in our problem. Therefore, the sliding 
mode surface of Equation 2.13 will be: 


=é+ke (2.14) 
= (hip — 41) +k(hiy — £1) 


To simpliy the computation, we use the lie derivative notion as 
follows: 


€=hires —hy (2.15) 
=hires — oN 
= hnnves — 2 (a) (2.16) 
Where ah, 


L h(x) = =" f(e) (2.17) 
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In order to find the sliding mode control law u, Equation 2.7 will 
be used: 


d 

FAG + ke) = — kpsgn(S) 
é+ké=—k,sgn(S) (2.18) 

(hip = Le hy = LgLrhyu) + k(hip == Lyhy) = kpsgn(S) 


Where 
Leh =2 
L4hy = — 2 +e(1—a7) xo (2.19) 
LgLshy =1 


Substitute Equations 2.19 onto equation 2.18: 


Ce — fo) + k(a1, — wa) = — kysgn(S) 
(hip — (—2, + (1 — 27)z2 + u))4 k(hip £1) = — kpsgn(S) 
(2.2 


Thus, the sliding mode control law u is obtained from Equation 
2.20 by: 


u= hay tar — €(1— 22 )22 + k(hip — 41) + kpsgn(S) (2.21) 


2.2.1.2 Solution Using MATLAB® 


The Nonlinear systems should be written in the following 
form—Sliding Mode Controller for a class of Nonlinear 
systems—State space equations 
x=f(x)+gu 


Number of states:=2 
Number of inputs:=1 
Parameters e 


parameters = 

e 

The vector f(x) :=[x2;—xlt+ex (1—-x17*2) *x2] 
The vector g(x) :=[0;1] 


The output variables:=x1l 
The sliding mode control law for SISO systems=: 
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Uc = 

d3yr + xl + kpxsgnS + klx«(d2yr — x2) + e*x2x*(x1°2 — 1) 
The Sliding mode surface:= 

Surf = 

d2yr — x2 + kl«(dlyr — x1) 

The derivative of Sliding mode surface:= 

dSurf = 

d3yr + xl + kl«(d2yr — x2) + exx2*(x1*2 — 1) 

2.2.2 DC motor angular position control 


The mathematical model of the DC motor [16] is given by the 
following state space representation: 


Ly, =22 (2.22) 
: b km 
t= zr + 7 08 (2.23) 
ke R 1 
eT ee ee 
hy = (2.24) 


Where x, and x2 are the angular position and velocity of the shaft 
respectively, x3 is the current of the armature coil. u is the input 
voltage of the armature coil. The parameters are the resistance 
and inductance of the armature coil respectively, ke and k,,, are 
the speed and torque constants, b is the viscous friction and J is 
the moment of inertia. We seek to control the angular position 
using sliding mode control. 


2.2.2.1. Solution Using Hand 
The sliding mode surface is defined as: 
S=(—+k)""e (2.25) 


Where the degree relative r = 3 in this example; the error e = 
Aires — hi, Equation 2.25 will be: 


S=é+2ké+e (2.26) 
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The control law is given by differentiating the surface S: 
S= €42ké+ ke (2.27) 


Using the notion of lie derivative the errors in Equation 2.27 
will be: 


e= Mire . hy 
e= hiref ae hy = hiref a Lyh(x) 
é =hires — hi = hirer — L7h(2) (2.28) 


@ =hirep — hi = hires — L¢h(x) — Lg LFh(a)u 


After having defined the derivative of the sliding mode surface, then 
Equation 2.7 will be used to construct the sliding mode surface: 


€ + 2ké + k?é = —kysgn(S) (2.29) 
Substituting Equation 2.28 onto 2.29 yields: 
hires — L¢h(x) — LgLFh(x)ut 2k(hirep —L3h(x)) — (2.30) 
+k? (hires — Lph(a)) = —kpsgn(S) 
The control law u is then given by: 
LgLh(x)u = hiret — L¢h(w) + 2k(hirep — L¢h(x)) (2.31) 
P k? (hiref | Lyh(z)) + kpsgn(S) 


Where 
Lyh(x) =x 
Li n(x) = — ory + ae (2.32) 
L3h(x) = rue — Anis) | 7 (Fm (Raz + het) 
LgLeh(x) = + 


Substituting Equation 2.32 onto Equation 2.31, the sliding mode 
control law is: 


sis ke R 
U= hires + — £2 + + #3 


L L 
+ 1 (bay —e kim23) 1 (Ras + kex2) 
+2 re b km 
(hires — G0) + 5 —_») 
(2.33) 
+k? (Airep — 2) + kpsgn(S) 
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2.2.2.2 Solution Using MATLAB 


The Nonlinear systems should be written in the following 
form—Sliding Mode Controller for a class of Nonlinear 
systems—State space equations 
x=f(x)+gu 


Number of states:=3 
Number of inputs:=1 
Parameters b J km ke LR 


parameters = 


b J km ke LR 


The vector f(x) :=[x2;—(b/J) *x2+(km/J) *x3; 
—(ke/L) *x2—(R/L) *x3] 
The vector g(x) :=[0;0;1/L] 


The output variables:=xl1l 
The sliding mode control law for SISO systems=: 


d4yr + kpxsgnS + kl«(d3yr + (b*x2)/7T — (km*x3)/J) 
+ k2*«(d2yr — x2) + (km*((R*x3)/L + (ke*x2)/L))/Jd 
— (bx ((b*x2)/JT — (km*x3)/d))/Jd 

The Sliding mode surface:= 


Surf = 


d3yr + k2*(dlyr — x1) + ki«(d2yr — x2) + (b*x2)/d 
— (km*x3)/J 


The derivative of Sliding mode surface:= 
dSurf = 
d4yr + kl«(d3yr + (b*x2)/J — (km*x3)/J) + k2«(d2yr — x2) 


+ (km«* ((R*x3)/L + (ke*x2)/L) ) /J— (bx ( (b*x2) /T 
— (km*x3)/J))/d 
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2.2.3 Permanent Magnet Synchronous motor speed control 


The dynamical model of PMSM [15] describing the electrical and 
mechanical part is given by: 


Be 2 Ly be 1 
ty = i ane T2X3 Te Ud 
R Ly ® 1 
to = 2.34 
v9 im Be) PTs ne in Ug ( ) 
Dy Ly —La f 1 

3 = xv LyX ax T 

3 yo? Pp J a io 


Where x1 and x2 are the d and q axis stator currents respectively; 
x3 is the mechanical speed of the motor; ug and u, are the d 
axis and q axis stator voltages respectively; R and Lg = La are 
the winding resistance and inductance on the d and q axis. J is 
mechanical inertia of the motor; 7 is the electrical torque. 

The objective is to control the mechanical velocity +3 and the 


x1 current. ae | h = | U1 (2.35) 


X3 


2.2.3.1 Solution Using Hand 


We note that the relative degrees of the current and the mechanical 
velocity are r; = 1 and rg = 2, respectively. Therefore, the sliding 
mode surfaces are: 


d 

Sy = (a om Rye = €1 (2.36) 
d : 

So = (a IF k)"2-1e, =é€o+ keg 


Where the errors are ey = Aipep — hy and e€2 = horet — hg. Thus, 
the control law is then obtained by differentiating the sliding mode 
surfaces in Equation 2.37 and equating to (2.6), yields: 


€, = —kpisgn(S1) (2.37) 
€2 + kég = —kposgn(S2) 
Using the notion of lie derivative Equation 2.37 will be: 
hires — Lphy — Lg, hiua = —kpisgn(S1) (2.38) 
horep — Liha — Lg, Lphota — Lg, L photg + k(hares — Lh) 
= —kpesgn(S2) 
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Rearranging Equation 2.38 and solve for the control law ug and tg 


gives: 
Ud _ Lgilshy 0 —t 
Uq Lg, Lhe Lg, L pha 
. hire = Lyhy =F kpisgn(S1) | 
hares — Leho + k(haret — Lph2) + kpasgn(S2) 
Where: 


Lyilphi 0 _ is 0) 
Loilphe Lg Lhe a pat2(La—Ly p&+pxi(La—La) 


(2.39) 


(2.40) 


ILa JLa 
Lphy = — TZ % + py tats 
Lyg—-L 1 
Lsh2 =p a Pp aa ora» ns ca 
T+ fas — p®x — *%1x2(Lq — Lg) 
L3hy =(f = ”) (2.41) 
= (p® + px; (Lg — Lq))(Rxz + p® + pLgx) £2) 
JLa 
_ pto(Lq — Lq)(Ra — pLgten3) 
JLa 


2.2.3.2 Solution Using MATLAB 


The Nonlinear systems should be written in the following 
form—Sliding Mode Controller for MIMO 
Nonlinear systems—State space equations 
x=f(x)+gu 


Number of states:=3 

Number of inputs:=2 

Parameters R Ld Lq phi f J p Tl 

The vector f(x) :=[—(R/Ld) *xlt+px (Ld/Lq) *x2*«x3; 

—(R/Lq) *x2—px (Lq/Ld) *x1*x3—p*phi/Ld; 

p* (phi/J) *«x2—px ( (Lq—-Ld) /J) «x1*x2 
—(£/J5)*x3-T1/J] 

The vector g(x):=[(1/Ld) 0;0 (1/Lq);0 0] 

The output variables:=[x1;x3] 

The relative degree of hl 

equal:=1 
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The relative degree of h2 

equal:=2 

The sliding mode control law for MIMO systems 
Is Given by Uc=:inv(Lhg) « (S) 
The function S:= 


S= d2yrl + kpl*sgnS1l + (R*x1l)/Ld — (Ld*p*x2*x3)/Lq ; 


d2yr2 + kp2*sgnS2 + T1/J + k1l«*(d3yr2 + ((px*phi)/J 
+ (p*xxl« (Ld — Lq))/J)*((R*x2)/Lq + (p*phi) /Ld 


+ (La*p*x1l*x3)/Ld) — (f*(T1L/J + (f£*x3)/d — (p*phix*x2)/Jd 
— (pxxl*x2* (Ld — Lq))/J))/J + (p*xx2* (Ld — La) * ((R*x1) /Ld 
— (Ld*p*x2*x3)/Lq))/J) + (£*x3)/7J — (p*phix*x2)/J 

( 


p*exl*x2x* (Ld — Lq))/J 
The Matrix Lhg:= 
Lhg = 


[ 1/Ld, 0] 
[ (p*x2%* (Ld — Lq))/(J*Ld), ((p*phi)/J 
+ (p*xl* (Ld — Lq))/dJ)/Lq] 


The sliding mode control law for MIMO systems=: 
Uc =Ld« (d2yrl + kp*sgnS +(R*x1)/Ld — (Ld«p*x2*x3)/Lq); 


Lax (d2yr2 + kpxsgnS + T1/J + kix*(d3yr2 + ((p*phi)/J 
p*xl* (Ld — Lg))/d)«((R*x2)/Lq + (p*xphi)/Ld 
La*p*x1*x3)/Ld) — (f*(T1/J + (f*x3)/7 — (px*phi*x2)/J 


(Tx 
ef 
#. 
— (p*xl*x2* (Ld — Lq))/J))/T + (p*x2* (Ld — La) *« ((R*x1)/Ld 
( 
( 
( 
( 


Ld*p*x2*x3)/Lq))/J) + (£*x3)/d — (p*phi*x2)/Jd 
pxxl«*x2* (Ld — Lq))/J))/(p* (phi + Ld*xl — La*x1)) 
Lq*x2«* (Ld — Lq)*(d2yrl + kp*xsgnS + (R*x1)/Ld 
Ld«p*x2*x3)/Lq))/(phi + Ld*xl — La*x1) 


2.3. TUNING OF SLIDING MODE CONTROL 
PARAMETER 


The parameters of the sliding mode control law should be cho- 
sen in order to enhance the performance of the controlled system. 
However, those parameters are said to be chosen arbitrarily using 
trial-and-error method, this method may take a longer time to find 
approximately the exact values. 

We will use a simple method that leads to approximately the 
exact values by identification with a characteristic of first or second 
order systems with the Equation 2.7. Let’s suppose that we have a 
nonlinear system with a relative degree r = 2; therefore, Equation 
2.7 will be: 


é+ ké = —kpsgn(S) 
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since the sgn(s) = is equation will be: 
é+ké=— bya 
5+ i; 
é+ké = —ky- a is (2.42) 


We note that the |5| is a positive function written in terms of 
the error, which can be chosen as a small value |S| << 1. Thus, 
Equation 2.42 will be: 


k k 
é+(k+—")é+ se=0 (2.43) 
[S| [S| 


Equation 2.43 takes the form of a second order system: 
3? + 26ws +w? =0 (2.44) 


By identification to Equation 2.44 to 2.48, yields: 


k 
(k+ isp) = 2¢w 
k 
is] = (2.45) 


The variable w and ¢ are chosen to meet specific requirements 


2.4 MATLAB-BASED FUNCTION: SLIDING MODE 
CONTROL 


MATLAB programs listed on this chapter are used to find the slid- 
ing mode control law for SISO and MIMO systems using Symbolic 
MATLAB Library. To simplify the use for these programs a test 
function is developed. After calling this function, the user has to 
enter the necessary functions which are: f(a), (g) and the out- 
put vector h(x) using symbolic variables. The principle MATLAB 
functions are: 


NonContSidFed: This function calculates the lie derivatives of 
the given SISO dynamical nonlinear system 


SlidingModeTerms: This function calculates the sliding mode 
terms like the surface and its derivative, as well as the control 
law U 
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The derivatives of the output y for SISO nonlinear systems are 
described by the following variables: 


ly yo ow. yt ]=[ dlyr dayr ... d(n+ yr ] (2.46) 


MIMOSlidingModeLieDer: This function calculates the lie 
derivatives for MIMO dynamical nonlinear systems 


MIMOSlidingModeController: This function calculates the 
sliding mode terms and control law for MIMO nonlinear sys- 
tems 


The derivatives of the output y for MIMO nonlinear systems are 
described by the following variables: 


Yi Yy we yt dlyr1 d2yrl .... d(n+1)yr1 

Yo Yo um» yf? dlyr2 d2yr2 ...., d(n+1)yr2 

Yn Yn we y® dlyrn d2yrn .... d(n+1)yrn 
(2.47) 


2.5 MATLAB PROGRAMMES DESCRIPTION 


The developed programmes are used to find the sliding mode 
control law for any types of nonlinear systems (SISO and MIMO) 
using MATLAB symbolic library. Hence, the user will put three 
separate MATLAB codes on the same folder which are: 
MATLAB name functions of SISO systems 


e Testfunction 
e NonContSidFed 
e SlidingModeTerms 
MATLAB name functions of MIMO systems 
e Testfunction 
e MIMOSlidingModeLieDer 
e MIMOSlidingModeController 


After calling the test functions by users, they will be asked to 
introduce the system functions f(x), g(a) and the output vector 
h(a) as symbolic functions. Therefore, each program from those 
listed above will output useful functions which are: 
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(1) The functions MIMOSlidingModeLieDer and NonContSidFed 
will give us for variables which are: 


(a) The lie derivatives of h along the vector field f and g which 
are represented by: Lhf and the decoupling matrix Lhg 


(b) The vector of all lie derivatives along the vector field g 


(c) The relative degree r of the systems 


(2) The functions MIMOSlidingModeController and SlidingMod- 
eTerms will give us for variables which are: 


(a) The MIMOSlidingModeController will provide as with: 
le, der, Sur f,dSurf,Uc, S], where e is the errors vector, 
der: represent the derivatives vector, Surf: is the sliding 
mode surface; dSurf: represent the derivative of the sur- 
face; S is a vector contain the switching terms such that 
the controller Uc is equal to Uc = inu(Lhg) * S 


(b) [Surf,dSurf,dd, K,Uc] =SlidingModeTerms(Hc, L,r, 
Lhg); Surf is the surface of a SISO system; dSurf: The 
derivative of the surface; K: The controller parameter vec- 
tor; Uc: the control law for SISO systems 


2.6 MATLAB CODES FOR SISO SLIDING MODE 
2.6.1 MATLAB test function for SISO systems 


clear all;clce 

disp (' me 

disp('The Nonlinear systems should be written in the 
following form'); 


disp ('—Sliding Mode Controller for a class of Nonlinear 
systems—'); 

disp(' State space equations x=f (x) +gu rs 

disp (' ')i 


n=input('Number of states:='); 
nin=input ("Number of inputs:='); 
x=sym(zeros(1,n)); 
u=sym(zeros(1,nin)); 


for j=lin 
eval(sprintf('syms x%d',4)) 
x(:,j)=sprintf('x%d',j4); 
end 
for k=l:nin 
eval (sprintf('syms u%d',k)); 
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u(:,k)=sprintf('utd',k); 


end 

syms ul 

f=input ('The vector f£(x):=','s'); 
g=input ('The vector g(x):=','s'); 
He=input ('The output variables:=','s'); 


[Lhf, Lhg, dh, L,u, r]=NonContSidFed(f,g,Hc,x,ul); 
[Surf, dSurf,dd,K,Uc]=SlidingModeTerms (Hc,L,r,Lhg) ; 


2.6.2 MATLAB function for SISO systems 


function [Lhf,Lhg,dh,L,u,r]=NonContSidFed(f,g,Hc,x,v) 

The NonContSidFed MATLAB function finds the lie 
derivatives and the decoupling matrix that has to be used 
to find the Sliding mode controller, by providing this 
function with: 


f£ : The system function f (x) 

g : The system output function g(x) 

x : The state vector x=[x1,x2,....,xn] 

Vv : The new input vector 

He : The output of the system hn(x)=xn 

The program has to give the user the following functions 
Lhf : The lie derivative of the vector h(x) along 


f(x) Lhf=[L* {r—1l}hfn] 


AHP AP AP AP oP AP AP AP AP AP AIP AP AP AIP OP AIP OP 


Lhg The lie derivative of the vector field h 
along the function g; Lhg=[Lg1Lfhn] 
L : The lie derivatives L=[Lhfl L*{2}hfl....... 
L*{r—1}hf1] 
x=f (x)+gu 
L=[]; 
Lhf=[];Lhg=[]; 


if (nargin<4) 
error('Not enough input arguments!'); 
end 
nl=length (Hc) ; 
n=length (x) ; 
k=1;kk=1; 
while (k<n+1) 
dh=LieDer (Hc,x,n); 
Lg=dhxg; 
L=[L,dh«f]; 
if (Lg~=0) 
r=k; 
Lhf=[Lhf£;dh«f]; 
Lhg=[Lhg; Lg]; 
break; 
end 
Hc=dhxf; 
k=k+1; 
end 
u=inv (Lhg) *« (Lhf—v) ; 
function dh=LieDer (Hc,x,n) 
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dh=[]; 
for ii=l:n 


d=diff(Hc,x(ii)); 
dh=[dh,d]; 


end 
end 
r=k; 
end 


function [Surf,dSurf,dd,Kk,Uc]=SlidingModeTerms (Hc, L,r,Lhg) 


ae oP ol? 


That takes 
xp=f ( 
y=h ( 


AHP AP oP AP oP oP AIP AAP AIP AP OP A? 
BK 


k=[]; 


The SlidingModeTerms MATLAB function finds the sliding mode 
surface and its derivatives for SISO nonlinear systems 


the following form: 

x) +gu 

x) 

The output of the system Hc=x 

The lie derivatives L=[Lhfl L*2hfl....... 
L* {r—1}hf1] 

The relative degree of the system r=rl 
Sliding Mode surface for SISO systems 
Surf=(d/dt+k) *(r—l)e 

The derivative of the sliding surface 
Surf=d ((d/dt+k) * (r—1) e) /dt 

The controller parameter 
vector K=[1 kl k2,....kn] 


syms kp Uc U sgnS 


if (r==1) 


k=sym(ze 
for jj=1 


ros(1,xr—1)); 
2-1 


eval (sprintf('syms k%d',4jj)); 
k(:,jj)=sprintf ('ked', 433); 


end 

K=[1 kl]; 
end 
dd=sym(zeros 
for ii=l:rtl 

eval (spr 


ddi(:,i1) 
end 


e=[dd(1)—He] ; 


%s=[dd(r)—Hc 
s=[]; 

k=1; 

for kk=2:r 


(1,r)); 


intf('syms d%dyr',ii)); 
=sprintf ('dtdyr',ii); 


, 


li 


e=[e; (dd(kk)—L(k))]; 


k=k+1; 
end 
LL=fliplr (L) 


, 
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dl=fliplr (dd); 


for kk=1l:lir 
s=[s; (dl (kk)—LL (kk) )]; 
end 
S=flipud(e); 
Sp=s; 
Uc=Lhg*U; 
disp(['The sliding mode control law for SISO systems= 
Uc=Kx*Sp+kp*sgnS 


disp(['The Sliding mode surface:=']); 

Surf=K*S 

disp(['The derivative of Sliding mode surface:=']); 
dSurf=KxSp 


f')); 


2.7 MATLAB CODES FOR MIMO SLIDING MODE 


2.7.1 MATLAB test functions for MIMO systems 


clear all;clc 


disp (' 
disp('The Nonlinear systems should be written in the 
following form'); 


disp ('—Sliding Mode Controller for MIMO Nonlinear 
systems—'); 

disp(' State space equations x=f (x)+gu 

‘)i 

disp (' 


n=input('Number of states:="'); 
nin=input ("Number of inputs:='); 
x=sym(zeros(1,n)); 
u=sym(zeros(1,nin)); 


par=input ("Parameters ','s'); 
eval(sprintf('syms %s',par)); 
parameters=sprintf('%Ss',par); 


for j=lin 
eval(sprintf('syms x%d',4)) 
X(:,j)=sprintf('x%d',j); 

end 

for k=l:nin 
eval(sprintf('syms u%d',k)); 
u(:,k)=sprintf ('usd',k); 


end 

f=input('The vector f(x):=','s'); 
g=input ('The vector g(x):='); 

h=input ('The output variables:=','s"'); 


f=sym(f) ;h=sym(h) ;g=sym(g); 
[Lhf Lhg L r]=MIMOSlidingModeLieDer (f,g,h,x); 
[e,der, Surf, dSurf, Uc] =MIMOSlidingModeController(h,L,r 


,Lhg); 
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2.7.2 MATLAB functions for MIMO systems 


function [Lhf Lhg L r]=MIMOSlidingModeLieDer (fx,g,h,x) 

This programme finds a nonlinear controller which is the 

input output feedback linearization controller where the 

system should be given in symbolic format this function 

works well with SISO and MIMO systems which is written in 


AP oP ol? 


% compact form like dax=f (x) +gu; 

% y=h (x); 

SX : The state vector x=[xl,x2,...,xn] 

Sh The output vector h(x)=[x1;x2;...;xn] 

& Ex The f(x) the describes the system 

am? | The g.u vector of the output system 

Su The output vector 

% Lhf The lie derivative of the vector file 

% h along the function f(x) 

% Lhf=[L* {rl—1}fh1;L* {r2—1}fh2;...; 

& L* {rn—1}f£hn] 

% Lhg The lie derivative of the vector 

% field h along the function g 

% Lhg=[Lg1Lfhl, Lg2Lfhl,..... ,LgnLfhl; 
Pee ea ’ 

% Lg lithn;, boeZEian; «44.0.3 , LgnLfhn] 

ae The degree relative of the system 

& ca gl ee a a avons | 

SL : The lie derivative vector 

% EH [LREW WED osc cacy Lhfn; 

% L*2htl DL 2he2...o.b 2heiny 

% Lo (ri—1) hel LL” (42-1) hee os ccs 

% L* (rn—1)hfn 

] 

Ee 


f nargin <4 
error('Not enough input argument'); 


end 
k=1;L=[];kk1=0; 
Lhg=[];Lhf=[];vc=1; 


nb=length (h) ; 
r=zeros(1,nb); 

while k<length(h) +1 

h1=h (k) ; 

for i=l:nbt1l 

df=Lie_ Derivative (hl,x); this Lie derivative 
funetien 

solve for the g 


Ae dP ol? 


[lhf lhg]=solve_lie_der (df,fx,g); 
L=[L, lhf]; 
[n b]=size(lhg); 
for ii=l:in 
d=any(lhg(ii,:) “= 0); 
end 
if d==1; 
disp(['The relative degree of h',num2str(k)]); 
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disp(['equal:=',num2str(i)]); 
r(1,k)=1; 
break; 
else 
hl=lhf; 


end 

if i==nb+1 && d==0 
disp(['The system dose not admit an input output 
feedback linearization']); 
return; 

end 


end 

Lhg=[Lhg; lhg]; 

Lhf=[Lhf; lhf]; 

k=k+1; 

end 

function df=Lie_Derivative (h, x) 

The LieDerivative MATLAB function is used 

to find the jacobian vector of a given output 


ale PoP oP lO 


if nargin<2 & nargin== 

error('not enough input argument'); 
end 
df=[]; 
n=length (x) ; 
for ii=l:n 

XX=x (ii); 

dff=diff (h, xx); 

df=[df,dff]; 
end 
df; 
end 
function [lhf lhg]=solve_lie_der (df, fx,G) 
This equation solves the Lie derivativies that is 
multiplied by 
The — vectors Lhg=Lh*g*u where 


AP JP oP oP oP oP al? 


LH : The Lie derivativies of h(x) along the vector 
field fx 

Laps : The vector field f(x) that describes the system 
G : The vector field of the input 

LHg=[]; 

lhf=df«fx; 

[n,b]=size(G); 

for ii=l:b 


Loh=df*G(:,ii); 
LHg=[LHg, Lgh]; 
end 
lhg=LHg; 
end 


h : Is the output function h(x)=[x1l;x2;...;xn] 
x : The state vector x=[xl,x2,..... Fe 2c] 
df : The jacobian of h along x 
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end 


function [e,der, Surf, dSurf,Uc,S]=MIMOSlidingModeController 


(h,L,r, Lhg) 
% The mimo nonlinear systems described by the following 
% xp=f (x) +gu 
Se : Is the error vector 
& e=[el e2 ....enj;del de2 ....den;..... ;dnel 
% dNe.... dNen] 
% der : The derivative of the output vector 
3 der=[hl h2 ...hnj;dlhl dlh2....dNhn;...... ;dNh1 
% dNe....dNen] 
$ Surf The Sliding Mode Surface for MIMO systems 
Es Surf=(d/dttk) *(r—-l)e 
% dSurf The derivative of the sliding mode surface 
% dSurf=d ( (d/dt+k) * (r—1)e) /dt 
= Uc The Sliding mode control law 
Sh The output vector y=[h1;h2;....;hN] for MIMO 
% systems 
6 L : The lie derivative vector of the MIMO system 
% L=[Lhtl LRE2 25.3 deco Lhfn; 
% Lr 2htl Lo 2hnt2... 2... 2hen; 
% Behe: Bent] ss den L*rhfn] 
ey ae : The relative degree vector r=[rl,r2,r3,....,rn] 
% Lhg : The lie derivative of h along the vector field g 
% Lhg=[Lg1lLfhl, Lg2Lfhl,..... ,ugnLfhl; 

fetes ee ean é 

% LolLfthn, Lo2bihny +... , LgnLfhn] 
e=[];kk=1; 
rr=max(r); 
der=[];S=[];Surf=[];Ss=[]; 


nb=length (h) ;d=[];qq=1;q=1; 
if (r==1)& (nb==1) 

error('The output vector h should be of length >=2'); 
end 


if nargin <4 
error('Not enough input argument'); 
end 
LL=sym(zeros (rr+1,nb)); 
der=sym(zeros (rrt+l1,nb)); 


syms Uc U 
LL(1,:)=h; 
k=[]; 


K=sym(zeros(nb,rr)); 
KK=sym(zeros (nb, rr+1)); 
KK(:,1)=sym('1'); 
shit, 1) =1; 
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for ll=l:nb 
R=r(1l); 
k=sym(zeros(1,R—1)); 
for jj=1:R-1 
eval (sprintf('syms k%d',4}4)); 
k(:,jj)=sprintf ('k%d', 44); 
end 
K(11,:)=[1 k]; 
end 
KK (:,nb:rr+1) =K; 
F=subs (KK(:,nbtl:rr+1),{1},{0}); 
KK(:,nb+1:rr+1) =F; 


%% The functions sgnSl, Sgn(S2),....sgn(Sn) 
% The parameters kpl, kp2,....kp3 


sgns=sym(zeros(1,nb)); 
kp=sym(zeros(1,nb)); 
for jj=l:inb 
eval (sprintf('syms sgnS%d',j4)); 
sgns(:,jj)=sprintf('sgnS%d', 44); 
end 
for jj=linb 
eval (sprintf('syms kp%d',jj)); 
kp (:, 33) =sprintf('kp%d', 33); 
end 


ole 
ole 


for jj=l:nb 
R=r (jj); 
dd=sym(zeros(1,R)); 
for ii=1:Rt1 
eval (sprintf('syms d%dyr%d',ii,kk)); 
dd(:,ii)=sprintf('d%dyr%d',ii,kk); 
end 
d=[d,dd]; 


kk=kk+1; 
end 


for jj=l:nb 
R=r(jj)3 

for mm=2:Rt1 
LL (mm, 33) =L (qq) ; 
qq=qqtl; 

end 

for nn=1:R+t1 
der (nn, jj)=d(q); 
a=qtl; 

end 

end 
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e=flipud(der—LL) ; 


for bb=1:nb 
Kd=fliplr (KK (bb,nb:rr+1)); 
ee=e(l:irr,bb); 
nm=kp (bb) *sgns (bb) ; 
sc=Kdxeetkp (bb) «sgns (bb) ; 
S=[S;sc]; 
Ss=[Ss;Kdsee]; 

end 


es=flipud(e); 

for cc=l:nb 
R=r (cc); 
dK=fliplr(KK(cc,nb:Rt1)); 
esf=(es(1:R,cc)); 
surface=dKx (esf); 
Surf=[Surf; surface]; 


end 
dSurf=Ss; 
Surf=Surf; 


[nh bh]=size(Lhg); 
disp(['The sliding mode control law for MIMO systems']); 


disp ([‘———Is. given by Uc=:iny(Lhg) « (3)———"]});; 
disp (['———The function $:=———————"]) ; 

S 

disp (['~———The Matrix Lhg:=———————']); 

Lhg 


for ii=1l:nh 
invg=any (Lhg(ii,:) ~=0); 
end 
if (invg~=0) 
disp(['The sliding mode control law for MIMO systems=:']); 
Uc=inv (Lhg) « (S) 


else 
disp(['The sliding mode controller is not possible']); 
disp(['for this output vector try to choose another one!']); 


end 
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